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The maximum efficiency in organic light-emitting diodes (OLEDs) depends on the ratio, r = 
ks/kr, where ks (kr) is the singlet (triplet) exciton formation rate. Several recent experiments 
found that r increases with increasing oligomer length from a value r w f in monomers and short 
oligomers. Here, we model exciton formation as a multi-phonon emission process. Our model is 
based on two assertions: (i) More phonons are emitted in triplet formation than in singlet formation, 
(ii) The Huang-Rhys parameter for this phonon emission is smaller in long oligomers than in short 
ones. We justify these assertions based on recent experimental and theoretical data. 

PACS numbers: 73.50. Jt,73.50.Gr,78.60.Fi 
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I. INTRODUCTION 

The maximum possible internal quantum efficiency, 
Vmax, of fluorescent-based organic light emitting diodes 
(OLEDs) occurs when the probability that the injected 
carriers form excitons and the quantum yield for singlet 
emission are both unity. rj max is then determined by 
(and identical to) the fraction, f s of injected electrons 
and holes (or negative and positive polarons, respec- 
tively) that pair to form emissive spin-singlet excitons, 
rather than noncmissive triplet excitons. If the process 
by which these excitons form were spin independent, then 
Vmax would be limited to 25% based on spin-degeneracy. 
However, recent reports indicate that Vmax in OLEDs 
ranges between 22% to 83% 0, H E @ . The 
exact value of rj max and the reason for this variation, 
however, have remained controversial. Indeed, even the 
notion that r) max can be larger than 25 % is currently not 
universally accepted 0. 



A. Overview of experimental results 

Two entirely different experimental approaches have 
been employed to study spin-dependent exciton forma- 
tion for OLEDs and thin films: 

(i) Experiments 0, that determine the 
singlet generation fraction fs directly in live OLEDs. 
For fluorescent devices typically only the singlet emis- 
sion can be measured, information on triplet density is 
missing and rather involved models have to be employed 
to obtain fs 0. Wilson et al. however, have recently 
shown 0] that in OLEDs made from organic semicon- 
ductors that exhibit spin-orbit coupling, the strong in- 
tersystem crossing implies that both singlet and triplet 
emission (fluorescence and phosphorescence) can be si- 
multaneously observed. This could be used to reliably 
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FIG. 1: Models used to describe polaron levels and opti- 
cal transitions depicted here for the positive polaron. (a) 
Electron-phonon (SSH) model, (b) Molecular orbital picture. 



determine fs by comparing the relative intensities of fluo- 
rescence to phosphorescence for optical excitation (where 
initially only singlet excitons are formed) with that for 
electrical excitation (where both singlet and triplet ex- 
citons are formed). Importantly, they found fs — 57 % 
for devices made from a Platinum-containing polymer, 
but fs — 22 % for the corresponding monomer OLEDs. 
This suggests that exciton formation is spin-independent 
for the monomer, but that a spin-dependent formation 
process is effective in the polymer. 

(ii) Experiments 0, IE E] that measure the ratio, 
r = ks/kx of the spin-dependent exciton formation rates 
for singlet and triplet excitons, respectively. Such experi- 
ments manipulate the spin state (using electron spin res- 
onance techniques) of the pairing polarons, and measure 
the effect on exciton formation rates. These experiments 
consider photogenerated polarons in the film and use the 
fact that antiparallel spin polaron pairs can either form 
singlet or triplet excitons, whereas parallel spin pairs can 
only form triplets. These optically detected magnetic 
resonance (ODMR) techniques are modulation experi- 
ments where the resonant /x-wave field is periodically 
turned on and off. Since the experiment is performed 
at low temperature, spin-alignment is conserved during 
the half- wave with /x-wave field off, and polaron recombi- 
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FIG. 2: Magnetic-resonance experimental data for the ratio 
r _1 —kr I 'ks of spin-dependent exciton formation rates in var- 
ious polymers and oligomers as a function of the peak photon 
energy of the Pi transition (lower x-axis). r _1 is also shown 
as a function of the inverse conjugation length 1/CL (upper x- 
axis) , which was determined from Pi (see text for discussion) . 
The line through the data points is a linear fit. *The Pi band 
of this polymer does not show a clear peak in the PA spec- 
trum, the Pi band extends to the longest wavelengths mea- 
sured. **The length of this oligomer was calculated. In addi- 
tion to the chemical names defined in the text, 3PE stands for 
the PPE trimer, PPE for poly(phenylene-ethynylene), Si-PT 
for silicon bridged polythiophene. For details consult original 
publications. 



nation/exciton formation obeys spin-statistics. However, 
during the half- wave with /z-wave field on, spin- 1/2 res- 
onance leads to rapid spin-flips of the recombining po- 
larons. Spin alignment is therefore not conserved, and 
each pair may choose whether to form singlet or triplet 
exciton. It can easily be shown 0, [|| that this leads to 
enhanced formation of the exciton with larger formation 
rate (leading to positive ODMR signal), at the expense 
of the more slowly forming exciton (that gives negative 
ODMR). In addition, the overall polaron recombination 
rate is enhanced, since the fast channel becomes allowed 
for all polaron pairs. Therefore changes occur in the pho- 
toinduced absorption (PA) from the triplet state, as well 
as the fluorescence from the singlet state upon magnetic 
resonance. In particular, from the /x-wave induced change 
in PA of the polaron pairs, r — ks/kx could be deter- 
mined 
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FIG. 3: The peak photon energies of the Pi polaron transi- 
tion in a variety of oligomers, namely solutions of (unsub- 
stituted) oligophenyls (OP, A, radical anion (RA)), alkyl- 
substituted (AS) oligophenylene-vinylenes (OPV, o, radical 
cation (RC)), alkoxy-substituted OPV (ffi, RC), end-capped 
oligothiophenes (OT, □, RC), films of AS OT (B, PA), AS 
oligothienylene-vinylenes (OTV, V, RC). The solid line is a 
fit to the data excluding ©. 



B. Polarons in 7r-conjugated semiconductors 

It is well-known that chemical doping or electrical 
charge injection results in the formation of polarons in n- 
conjugated semiconductors El- Fig- fflshows a compari- 
son between different models that have been used for de- 
scribing polarons in 7r-conjugated semiconductors. Panel 
(a) depicts the electron-phonon (e-p) or Su-Schrieffer- 
Heeger (SSH) model E3, EL E2 ■ It predicts that the e-p 
coupling causes a gap between valence and conduction 
band. In the singly charged system two localized polaron 
levels appear inside the gap. Experimentally one finds 
two optical transitions |13| that are interpreted as the 
Pi and P2 transitions. Panel (b) depicts the molecular 
orbital picture where HOMO and LUMO are the high- 
est occupied and lowest unoccupied molecular orbitals, 
respectively. 

The most generic model for polarons is the molecular 
crystal or Holstein model El ! it i s however not expected 
to be applicable in a quantitative way to 7r-conjugated 
polymers. This model yields for the "Pi" transition [Tl| : 
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Here A quantifies the e-p coupling strength (e.g. in 
eV/A),Mis the ionic mass and u>e is the Einstein phonon 
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frequency, W is the band width before inclusion of e-p 
coupling. The term in the bracket is the energy, V as- 
sociated with the e-p coupling 0|. Furthermore, the 
following two equations hold for the polaron binding en- 
ergy, Eb polaron j 

and the deformation or relaxation energy 

Er elax , polar on j 

respectively. 



E 



b,polaron ~ ^ ^1 



E, 



r elax , polar on ^ 1 



(2) 
(3) 



Il-conjugated oligomers are often used as model com- 
pounds instead of 7r-conjugated polymers because they 
can be obtained with a well-defined chemical structure. 
Although the molecular weight of polymers is typically 
much larger than that of oligomers, nevertheless it is es- 
tablished that the polymer should be viewed as a string 
of effectively independent segments, separated by chem- 
ical or physical defects. The length of these segments 
is called the conjugation- length (CL). Several recent cal- 
culations found that the energy associated with the e-p 
coupling decreases with increasing oligomer size. In par- 
ticular, Devos and Lannoo found that V = const /N in 
a more or less universal manner in acenes and fullerenes 
with various numbers, N of 7r-bonds and Shuai et al. 
im found that Eb jPO iaron in oligophenyls (OP) decreases 
with the number of phenyl rings in a similar manner. We 
may therefore state that, in average, that effect of e-p 
coupling significantly decreases with increasing oligomer 
length. 



C. Results of magnetic resonance experiments 

Using ODMR, it was found that r is a monotonously 
increasing function of the conjugation- length (CL), and, 
by extrapolation, that r ~ 1 for small molecules and 
monomers |(| . Electroluminescence and magnetic res- 
onance experiment therefore lead to the same qualita- 
tive conclusions, namely that exciton formation is spin- 
independent for the monomer, but that a spin-dependent 
process is effective in the polymer. Fig. [5] shows the ra- 
tio r =kr j kg of spin-dependent exciton formation rates 
in various polymers and oligomers as a function of the 
peak photon energy of the P\ transition (lower x-axis) 
obtained by ODMR spectroscopy |(| . In the original pub- 
lication the experimental dependence r(P{) was however 
reinterpreted in terms of the CL of the polymer films 
studied. This was possible, because it is known that 
there exists a (almost) universal relationship between P\ 
and the material's CL (see Fig. r _1 is therefore also 
shown as a function of 1/CL (upper x-axis), which was 
determined from Pi . 
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FIG. 4: Charge-recombination electronic couplings into sin- 
glet and triplet excited states in 6PV. Panel a) cofacial ar- 
rangement, panel b) head-to-tail arrangement. The charge- 
transfer state occurs at « 3.7eV. The data points were taken 
from Ref. 



D. Recent theoretical results 

Primary excitations in these materials are generally be- 
lieved to be excitons with a binding energy in excess of 
kT, where k is the Boltzmann constant and T the tem- 
perature. As a result of electron-phonon and electron- 
electron interactions, the lowest singlet (Si) and triplet 
(Ti ) excitons posses both different energies (the Si — Ti 
energy difference; the exchange energy K for the lowest 
excitations has been measured 0] to be K — 0.7eV in 
a variety of conjugated polymers) and different spatial 
wavefunctions (with T± displaying a more spatially con- 
fined character). 

Our model is based on recent work of Beljonne et al. 
[l7|. They developed a theoretical model to describe in- 
termolecular charge recombination in conjugated mate- 
rials. In their treatment, Beljonne et al. found it nec- 
essary to consider two configurations of the two polymer 
chains involved in the exciton formation/polaron recom- 
bination process, namely cofacial and head-to-tail (see 
Fig. They found that in the cofacial arrangement, 
by far the largest matrix elements are calculated for the 
lowest singlet Si and triplet T\ excited states (see Fig. 0] 
a). Similar results were obtained also by Tandon et al. 
jl9j . The situation is quite different for the head-to-tail 
configurations, where a number of different singlet and 
triplet excited states show significant electronic couplings 
to the polaron pair states (see Fig. ^b). The following 
picture has therefore emerged: Based on the electronic 
coupling, polaron recombination is a direct transition 
predominantly to the lower lying exciton states. Such 
transitions, however, have to pay a high price, since the 
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multi-phonon emission necessary to conserve energy has 
very low probability for reasonable values of the Huang- 
Rhys parameter. Therefore, unless the two chains are 
in exact cofacial arrangement, the exciton states formed 
with highest probability may not be the lowest exciton 
states. We may therefore adopt a picture where exciton 
formation has highest probability for an intermediate ex- 
citon state. 



the exciton formation rate has maximum probability av- 
eraged over the ensemble of configurations. This step is 
justified in more detail in Appendix A. With this defini- 
tion in mind we may write: 



kcr.S F Qvs 
r = — = rw 



"CR,T 



F 



(8) 



E. Theoretical approach 

In the case of charge-recombination (CR) processes, 
the semiclassical expression for the CR rate writes, within 
the Franck-Condon approximation, as |17|: 



^ / (AG + A s + vfrujph 

x 2_^F 0v exp - — 



1/2 



(4) 



4A 5 fcT J 

Here \ipi > and \tpf > are the wavefunctions of the ini- 
tial and final states, respectively, W is the perturbation, 
huiph is the energy of the most strongly coupled (optical) 
phonon, Fq v is the Franck-Condon factor of the tran- 
sition with zero and v phonons in the initial and final 
states, respectively, G° is the difference in free energy 
between intial and final state, and As is the (external) 
reorganization energy. Next we approximate the above 
expression, since the largest contribution will result from 
transitions that conserve energy, namely those for which 
AG + As + vfuo p h » 0. This can always be achieved, at 
least approximately, through emission of a number ve of 
phonons: 



where the subscript S and T denote the "effective" sin- 
glet and triplet states, and vg and vt denote the number 
of phonons required for energy conservation to form the 
S and T states, respectively. The definition of rw fol- 
lows by comparison of Eq. [S] to Eq. [7| and is essentially 
the ratio of the electronic matrix elements for singlet and 
triplet formation, respectively. 

Next we discuss the calculation of F „, where we 
closely follow the treatment by Barford et al. [2(j 
Inter-molecular interconversion is an iso-energetic pro- 
cess which occurs from the lowest vibrational levels of 
the initial polaron pair state to the final, intra-molccular 
exciton states at the same energy as the initial level. In 
actuality, the exciton formation process involves two con- 
formational transitions, namely the transition from the 
polaron lattice conformation to the exciton conformation 
in chain 1 and from polaron to groundstate conformation 
in chain 2. We therefore need to generalize our treatment 
in the following way: 
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We therefore obtain for kcR- 
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The ratio, r = |^ is therefore given as: 



J2 ns k CR,S ns 

J2n T k CR,T nT 



(5) 



(6) 



(7) 



where the sums extend over all singlet states and triplet 
states, respectively. 

For the sake of simplicity, we now replace the sums 
in Eq. [3 by a single "effective" state S and T, respec- 
tively. This effective state may be loosely identified with 
the intermediate exciton level discussed above, for which 



where F^) and Fq„ 2 are the Franck-Condon factors 
associated with the vibrational wavefunction overlaps of 
chains 1 and 2, respectively. 

Next we want to write out the expression for the 
Franck-Condon overlaps using the displaced oscillator ex- 
pression Fo„ = (e~ s S u } jv\. A useful simplification to 
this expression arises by noting that the geometric dis- 
tortions of the polarons and exciton polarons (namely 
the \ Y B U or 1 3 B U states) from the ground state structure 
are very similar. [2l| Thus, the Huang-Rhys parameter, 
Si for the l 1 ^ and 1 3 B U states relative to the positive 
polaron is quite small. Therefore, to the lowest approx- 
imation, only vi = has a non-negligible contribution, 
and S total ~ S2 = Sp, where the P stands for polaron 
conformation. In the next higher approximation, we may 
also include the term v\ = 1. As we show in Appendix 
B, if Si -C 1 and S± <C Sp, then we may combine the 
two Franck-Condon factors into a single one with Huang- 
Rhys parameter, S = Si + Sp. 

Therefore: 
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FIG. 5: Experimental values for Pi polaron transition energy 
in a variety of oligomers, namely solutions of (unsubstituted) 
oligophenyls (OP, radical anion (RA)), alkyl-substituted (AS) 
oligophenylene-vinylenes (OPV, radical cation (RC)), end- 
capped oligothiophenes (OT, RC), films of AS OT (measured 
by photoinduced absorption), AS oligothienylene-vinylenes 
(OTV, RC) together with calculated values for the polaron 
relaxation energy E re i aXtPO i aron in OPV. The experimental 
data for Pi were taken from Ref, whereas the calculated val- 
ues for Eb, P oiaron were taken from Ref. 
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Eas.llllandll2lare the final result of our model. 



II. COMPARISON BETWEEN EXPERIMENTAL 
RESULTS AND MODEL 



and E re i aXtPO i aron in 7r-conjugated polymers. However, 
the Holstein model is formulated for the infinite system, 
and does not consider finite size effects. It is less than ob- 
vious that Pi oc E re i aXtPO i aron holds also in oligomers. In 
particular, since we are interested in phonon emission, we 
need to distill finite size effects on the e-p coupling energy 
from finite size effects on the electronic energies (quan- 
tum confinement energy). Finite size effects on the e-p 
coupling energy have been studied by Devos and Lannoo, 
who found that V = const /N in a more or less universal 
manner in acenes and fullerenes with various numbers, N 
of 7r-bonds 01- We therefore expect that E re i ax ^ po i aron 
decreases strongly as the oligomer size increases. 



A. The conjugation-length dependence of the 
polaron relaxation energy 

To the best of our knowledge a direct measurement 
of E re i ax polaron has not yet been performed in pi- 
conjugated polymers and oligomers. However, Bredas 
and coworkers have calculated 15] E re i aXtPO i aron for 
phenyl-capped PPV oligomers as a function of the num- 
ber of rings N against which we can compare our exper- 
imental results for P\. Their results are shown in Fig. El 
Fig. [S] also shows a plot of the Pi polaron transition en- 
ergy in a variety of oligomers. It is seen in Fig. [5] that 
both experimental data for P\ and theoretical calculation 
of E re i ax polaron follow a very similar dependence on CL, 
as a matter of fact Pi oc E re i aXiPO i aron to a high degree of 
accuracy. This is strong evidence for the notion that Pi is 
indeed a measure of E re i ax<po i aron . Therefore both exper- 
imental and theoretical data find the that E re i aXtPO i aron 
monotonously decreases with increasing CL. In particu- 
lar (see Fig. [3J) , 
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+ k x 



1 

CL 



(13) 



Through comparison of the data shown in Figs. [21 and 
0it is therefore evident that the experimental results can 
be brought in agreement with Eci. llll onlv if Vj< — vg = 1, 
therefore 



Guided by the magnetic resonance experimental data 
shown in Fig. [21 we now study the CL dependence of 
Eq-El The CL-dependence of rw was studied previously 
by Beljonne et al. ^3 They found that rw depends only 
weakly on the CL, in particular this dependence is not 
strong enough to account for the material dependence 
of r shown in Fig. [2 We must therefore study the CL 
dependence of the phonon emission term in Eq. 1111 

There is an obvious connection between our model and 
the experimental magnetic resonance results. These re- 
sults (see Fig. |2J) are given as a function of Pi, whereas 



S'i 



r w v T 



(14) 
(15) 



Importantly, since the experimental relationship 
Pi (CP) is universal, then it follows naturally that r(CL) 
is universal. 



B. The conjugation-length dependence of the 
optical Huang-Rhys parameter 



Holstein model relates the two quantities through Eq. [3 
A similar relationship should therefore hold between Pi 



Our model is based on the assertion that the Huang- 
Rhys parameter is substantially smaller in polymers than 
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FIG. 6: Experimental values for the optical Huang-Rhys pa- 
rameter, S^g , in a series of OPV and OT. 



oligomers. Thus far, we have presented three justifica- 
tions for this statement: 

1. The optical polaron transition energy Pi strongly 
scales with CL, specifically as Pi = P\ t00 +kp 1 xCL. 
Furthermore this scaling is universal. 

2. A very similar scaling law was found by calcula- 
tions of E re i aX}PO i aron vs. CL in OPV performed 
by Bredas and coworkers. 
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FIG. 7: A fit of the experimental data for r" 1 using our multi- 
phonon emission model and together with the experimentally 
determined optical Huang- Rhys parameter, S^^. 

experimentally determined optical Huang-Rhys parame- 
ter, S±ib ■ Specifically we used the fit function 



r w v T 

It is seen that excellent agreement is achieved. 



3. The electron-phonon energy, V was found to scale 
as V oc 1/-/V, where N is the number of 7r-atoms 
in a large class of conjugated compounds, largely 
independent of the specific chemical structure. 

Here we now want to add another set of data to justify 
the main assertion on which our model is built. The lat- 
tice relaxation of the polaron state is expected to be sim- 
ilar to that of the lowest singlet exciton state l l B u . The 
corresponding Huang-Rhys parameter, S\ib u can be eas- 
ily measured by optical absorption and emission. We may 
therefore test our assertion by determining S\b u from ab- 
sorption and emission spectra of oligomers. 

Fig[f)|shows experimental values for the optical Huang- 
Rhys parameter, Sx^b u i m a series of OPV and OT. 
The exp erimental data were taken from Ref. [2^] and 
Ref. [23 • It is seen that S\ib u decreases with increasing 
CL, namely S\i b u = Si 1 b u .oo + k x tjj; m agreement with 
our expectation based on E re i aXtPO i aron . 

Fig. [7| shows the fit of the experimental data for 
using our multi-phonon emission model together with the 



III. SUMMARY 

We developed a model of spin-dependent exciton for- 
mation in OLEDs. We calculated r = kgfkr based on 
a multi-phonon emission process. The resulting equation 
for r therefore strongly depends on the Franck-Condon 
overlap integrals that are parameterized by the Huang- 
Rhys parameter. Guided by recent magnetic resonance 
experiments, we studied the dependence of Huang-Rhys 
parameter on the conjugation-length. We used two dif- 
ferent approaches: 

1. We relate the Huang- Rhys parameter to the 
polaron relaxation energy. Then we relate 
Ereiax, polaron to the optical polaron transition Pi, 
which has been measured experimentally in many 
oligomers. This procedure leads us to conclude that 
S = S 00 ~\~ k x -j=r£ . 

2. We relate the Huang-Rhys parameter in our model 
to the Huang-Rhys parameter that can be mea- 
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sured by photoluminescence spectroscopy. Such 
data are available in OPV and OT. This procedure 
also leads us to conclude that S = Soo + k x . 

Our work therefore leads to the following picture of 
exciton formation: Since the triplet exciton states lie 
lower in energy than singlets, more phonons must be 
omitted (required by energy conservation) for triplet for- 
mation than singlet formation. Since polymers have a 
small Huang-Rhys factor, then the emission of many 
phonons is unlikely, thus favoring singlet formation. In 
short oligomers, however, the Huang-Rhys factor is quite 
large, phonons are emitted easily, and singlet and triplet 
formation both become likely. 



IV. APPENDIX A 

The definition of a single effective state is certainly 
appropriate in the cofacial configuration, since only the 
lowest exciton states contribute to the exciton forma- 
tion process. In the head-to-tail configuration, however a 
whole series of states contribute. Fig.0]b) shows that the 
electronic couplings are roughly equal in all states that 
contribute, such that we may write: 



Therefore the effective state in the head-to- 
tail arrangement has an electronic matrix-element 
(ipi}W\ij}f :Tl ), and Franck-Condon Factor equal to 
^2 n Fov E , n > an d nas a binding energy VEtuJph, such that 



V. APPENDIX B 



Including only the v\ — and v\ = 1 terms we then 
obtain 



Fqv e — 



P (l)p(2) , p(l)p(2) 
a -S lo -S_ 



Sp QVE — l 
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We - 1)! 



-(Si+Sp) 



(S P B + v E S l Sp s ) 



v E \ 



(18) 
(19) 

(20) 



On the other hand, the Franck-Condon factor, Fo Vmax 
for a (hypothetical) Huang-Rhys parameter, S = Si + Sp 
equals 



52^m\w\ii> f , n )\*' 1 



4n\ s kT 



1/2 



Fo» E , n (17) 



2n 



T m\w\M\ UnXskT 



i 



1/2 



0v E . 



-(Si+Sp) 



(Si + Sp) 1 



F 0VB - - ( 21 ) 

We see that if S\ -C Sp, the two expressions are equal 
to first order. 
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